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We re-derive a formula relating the areal and luminosity distances, entirely in the framework of 
the classical Maxwell theory, assuming a geometric-optics type condition. 



I. INTRODUCTION 

It is accepted since Robertson [1] that the luminosity distance relates to the areal distance [2] as follows, in 
Friedman-Lemaitre-Robertson- Walker (FLWR henceforth) spacetimes, 

D=(l + z)R. (1) 
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Formula (1) is usually derived in a way that mixes quantum mechanical and classical concepts [3]. Robertson has 
shown this result within the framework of the classical Maxwell theory using conservation laws. His approach suffers 
. | from a number of minor erronous statements and gaps in the reasoning. One of that concerns the form of the 
i electromagnetic energy-momentum tensor, another the energy flux. A comment on the latter. The Maxwell equations 
in FLRW spacetimes possess a mathematical strictly conserved energy which is different - as it is quite common in 
PsJ ■ curved spacetimes - from the 'proper energy', the total energy detected by a comoving observer. The latter is directly 
related to the luminosity, but is not conserved; in fact it is equal to the former energy divided by the conformal 
factor (see below). Robertson assumes that the corresponding energy fluxes scale in the same way as the energies 
themselves; that can be true only approximately (but including all cases of astrophysical interest), as will be shown. 

These shortcomings do not influence the validity of the main conclusion of [1]. The aim of this paper is to provide 
a derivation for the cosmological distance formula which is consistently classical and fills some gaps in the proof of 
Robertson [1]. The order of the rest of this paper is following. Sec. 2 brings fundamental definitions and equations. 
Sec. 3 is dedicated to the derivation of an explicit solution of the Maxwell equations. Next section shows how one 
incorporates initial data into this closed form solution. The first part of Sec. 5 contains a proof that the electromagnetic 
energy momentum tensor can satisfy assumptions made by Robertson. Sec 5B accomplishes the proof of (1). 
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II. EQUATIONS 

The space-time geometry is defined by the FLWR line element, 

ds 2 = a(r?) 2 (-<V + dr 2 + u 2 dVL 2 ), (2) 

where 77 is a (conformal) time coordinate, r is a radial coordinate that coincides with the areal radius, dil 2 = 
dd 2 + sin 2 9d(j) 2 is the line element on the unit sphere with < < 2it and < 6 < n. 0(77) is a conformal factor and 
a is given by sinh(r) (for negative curvature k = — 1), r (for flat spacetime k = 0) and sin(r) (for positive curvature 
k = 1 ). Throughout this paper c and G, the velocity of light and the gravitational coupling constant respectively, 
are put equal to 1. The standard cosmological time is related to 77 by adr/ = dt. 
The Maxwell equations read 

V M F"" - 0, (3) 

where F^ v = d IJ- A u — d v A^ and is the electromagnetic potential. It is convenient to assume Aq = and the 
Coulomb gauge condition V,-A* = 0. In such a case there are two independent degrees of freedom, represented by the 
magnetic or electric modes. 

The magnetic modes are defined as follows 
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Z_1 m — 1 y/Hl + T) \ sm9 
-ism9d e Y lm d^j. (4) 

Here Yi rn are the spherical harmonics. The multipole expansion coefficients fi m depend on the conformal time r\ 
and the radial coordinate r. Straightforward caclulation shows that the equations of motion reduce to a system of 
hyperbolic equations 

(-d 2 + d 2 r )f lm= l l±±lf lm (5 ) 

for the multipoles fi m (v, r )- 
The electric potential read 



A = Sg^U-t ^777^1) yVW+ljhl m (v, r)Y lm dr + 
+h m {r),r)d e Yi m d6 + 

+ki m (r),r)imY lm d(j)j. (6) 

Here hi m (r],r) and ki m (r),r) arc multipoles. From the Coulomb gauge condition \7iA l = we have the relation 
between multipoles ki m (r],r) and hi m (r],r): 



+ l)h m = d r (h lm ) (7) 

where hi m = a 2 hi m . One can show that the first Maxwell equation (V^F^ = 0) is identity. The second one, 
(V ^F^ r = 0), reduces to the same hyperbolic equations as in Eq.(5): 

{-dl + d 2 r )h lm = l ^fih lm - (8) 

The two remaining Maxwell equations (V^F^ = where k = 6,<fi) yield 



a 

It's easy to show that (8) is equivalent to (9). In order to see this, differentiate both sides of (8) with respect r and 
use the Columb gauge condition (7); a quick calculation gives (9). 



III. COMPACT FORM OF SOLUTIONS 



Equations (5) and 8) have identical form; they can be written compactly as 

(-d 2 + d 2 r )cb lm = l -V±^<t> lm , (10) 

where <p = f, h. Below we will drop the irrelevant index m. The forthcoming constructions apply equally well to 
magnetic and electric degrees of freedom. The form of Eq. (10) suggests that each of the three cosmological models 
should be discussed separately, but we show later, that there exists a unique approach that works in all these cases. 
In this section we construct a solution of Eq.(10). 

Define a set of generating functions (as we will sometimes refer to them in the forthcoming text) of the form 

f(v,r)=f(r-n) (II) 

and also 



2 



g(v,r) = g(r + T)). 



(12) 



Lemma. For any / and g the function 



<f>i(r,v)=° i d r -d r ----dr(t±?-) (13) 
a a a 

V v ' 

I 

solves Eq. (10). 

Proof of Lemma. 

The proof uses the method of matematical induction. One can write the right hand side of (13) as 
a l d r | (§i=t) — ■ ■ ■ d r ( — - — -) j (herein we simply put after the first sign of differentiation). Thence one ob- 



l-i 

tains a recursive formula 



where 



d (7 

(r, 7]) = a r (J _i) (r, rj)-l— <^_i) (r, 77) (14) 



4> { i-i)(r,v) = o 1 - 1 d r -d r ----d r (ttl). (15) 

a a a 
. ' 

i-i 

For I = 1 we can write the right hand side of (13) as 

4>i — adr^—-^-- (16) 
a 

This satisfies our equation (10), as can be easily checked. Now let (j>i be a solution of: 

(-d, 2 +d r 2 )<p l = l -V±V<t> l . (17) 



We claim that <^n+i) defined as in (14) in terms of (j>i satisfies the equation (17) with I + 1 being put in place of /. 
That is easily proved by a direct calculation using the relation (14)), which shows that indeed {—d n 2 + <9 r 2 )0( (+1 ) = 

_ , . ' 

Let us remark that explicit solutions of Eq. (10) are known in the literature [4]. Eq. (10) has the same form as in 
the Minkowski spacetime (strictly saying, it is so when k = 0. If k = 1, —1, Eq. (10) can be cast into the Minkowskian 
form by a suitable change of variables (see for instance [5], [6] and [7]). The specific form of the solution, presented 
in Eq. (13), can be new. 

We would like to point out that the pair (/, g) should be specified independently for each multipole <f>i ; the related 
pair of functions will be labelled by the multipole number I. The initial data for the wave equation (10) are 0/, <f>i. In 
the next section we show that those initial data suffice to determine the pair (fi,gi), which in turn determine solutions 
of Eq.(10) in the whole spacetime. 

IV. GENERATING FUNCTIONS AND INITIAL DATA 

Below we show how, having initial data one can construct the generating functions / and g that appear in the 
formula (13). 

Let initial conditions be 

Mr,v)l\ v =o = Mr,0) d n Mr,v)i\ n =o = Mr,0) (18) 

Assume that a support of initial data is compact and contained within (a, b). Integrating l-times both sides of (13), 
one arrives at 

{9i + fi)\ v =o = Hi(r)+M l (r), (19) 
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where quantities Hi and Mi are linear combinations of multiple integrals 



Mi = M a(r) + M s a{r) I dna(n) 

rri rr e -i 

/ dr 2 <j{r 2 ) ■ J dr s a(r s ), 

J a J a 

Hi(r) = a(r) / dricr(ri) / dr 2 (j(r 2 ) 

J a J a 

'"Ja 



Here M's are constants of integration. 

From (13), taking into account (2) we have 



where 



rr rri 

h{r) = / dna{n) \ dr 2 (r{r 2 ) 

J a J a 

"Ja dna{ri) i ^W" r 

pr 

Ni = n + ^2n, d ri a(n) 

8=1 Ja 



rri rr 3 -i 

/ dr 2 a(r 2 ) ■ ■ ■ dr s a 

J a J a 



and iV's are the integration constants. The substraction of (22) from (19) yields 



(20) 



Mr, 0) = a 1 d r -d r - ■ ■ ■ d r dr{91 fl) . (21) 
a a a 

" v ' 

I 

Therefore performing an adequate number of integrations, one arrives at 

(/* - 9l)\ v = =-h-N u (22) 



(23) 



9i\ v =o = \{Hi{r) + Ii(r)) + ^(r) + ^N t (r) (24) 
Similary, adding (22) to (19) gives 

/i|„=o = l(Hl(r,0) J,(r,0)) + \m x \n x . (25) 

Only quantities Hi and // carry information depending directly on initial data. In contrast, let us stress out that 
the constants M's and TV's are free - hence M; and ./V; are not determined by initial data. Therefore the generating 
functions are not specified uniquely by initial data. On the other hand, inverting this argument one can say, that this 
nonuniqueness does not influence initial data. Since a solution <j>i of Eq. (10) is uniquely determined by its initial 
data, one infers that properties of the evolving waves do not depend on specific values of the constants. Therefore we 
are free to adopt such values that are convenient to us, and this fact justifies the choice that will be made at the end 
of this section. 

The quantities Mi and iVj are linear combinations of multiple integrals. Using the additivity of the integrals and 
splitting f a = + f b , one obtains 
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* 1 rr pn 

Mi = M a{r) + V M s a(r) / d ri a{n) / dr 2 
»=i • /6 



cr(r 2 ) 



dr s cr(r s ), 



No + y^Ns dna(n) dr 2 a{r 2 ) 

Jb Jb 



s=l 

■jf dr s a(r s ), (26) 



with some constants M's and TV's. In a similar token (but under additional condition r > b) one obtains 



* 1 pT pT\ 

Hi(r, 0) = h a + V]/i s cr(r) / dno-(ri) / dr 2 a(r 2 ) 
~s Jb Jb 



s=l 

y dr s a(r s ) 
I 



- rr rri 

Ii(r, 0) = i + V] ^ / dricr(ri) / dr 2 a(r 2 ) 
s=l Jb Jb 

••• / dr s( x(r s ). (27) 

Here /i's and i's are fixed constants, that depend on initial data, but whose specific form is not relevant. We would 
like to stress, that here the condition r > b is essential. 

We choose the integration constants M's, N's such that the right hand sides of Eqs (24) and (25) become nullity 
for r > b. Therefore the generating functions fi,gi vanish in the interval (6, oo) and become linear combinations 
(MaiAT )/2 in (0,a). 

V. THE LUMINOSITY FORMULA 

A lesson that can be drawn from the preceding section is this: an outgoing wave pulse that is initially comprised 
within (t\e — A, te) (as we assume from now on) will remain inside (t\e + r\ — A, te + vj). A similar statement can be 
formulated about ingoing waves. This provides an explicit proof of the Huygens principle [8] in Friedman cosmological 
models [9]. 

In what follows we will restrict our attention only to the outgoing pulse of radiation having a width A, and located 
initially within (r^ — A,r^). The stress-energy tensor of the electromagnetic field reads = F^F" 1 — (l/4 : )g^F 7 sF' yS 
and the time-like normal to a Cauchy hypersurface is (n p ) = (—a, 0,0,0). Define 

m(r, t) = - f dVaT° = 

JV(r) 

2ir jf drEj^^z) 2 + {d r 4if + ^j^ 2 )- (28) 

where dV is the proper volume element and the volume V(r) extends over the support of a radiation pulse that is 
enclosed inside the coordinate sphere S(r). t is the cosmic time, related to the conformal time by dt = adr]. One finds 
that 

47T 

d t m(r,t) = — H^iMidrfo- (29) 

The right hand side of Eq.(29) can be interpreted as the "energy flux" through the comoving coordinate sphere S(r). 
Assume that A << <j{te) (we will refer later to this condition as to the geometric optics approximation). Then it 
follows from expression (13) that one can approximate cf>i by its leading term, d r f ', roughly speaking (see Appendix 
for a sketch of the precise argument). Therefore one has approximate equality valid almost everywhere 
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dt m(r,t) = — Y.r = M +1 fif- (30) 
The total energy m = lim m(r, t) is strictly conserved, because the energy flux vanishes outside a coordinate sphere 

r — >oo 

S(ro + v(t))- There is an interesting observation that can be made about the outgoing electromagnetic radiation 
in the geometric optics limit. If one adopts our point of view in which the emitter is not a point but rather a big 
sphere (say, a surface or an envelope of a galaxy of a radius te — A), then one can show that |Tb r |,Tbo and T rr are 
much bigger that the remaining components of the energy-momentum tensor. In order to shorten the mathematical 
formulae, we will operate with quntities integrated over the coordinate sphere te- (In the pointwise treatment one 
would have to deal with expressions of the type ^ v Y/o^z'o^zi' , where the symmetric tensor Aw is given by Au> = 
(l/2)(d r (j)id v (j)i> +d r (f>i>d v (f>i) or A w = (l/2)(<9 r (^d r <^/ -d^i'd^i).) The dominant parts of j s d 2 S\T 0r \, J s d 2 ST 00 and 
J s d 2 ST rr behave like dr^id^i; this can be be approximated, in the limit of geometric optics, that is A << te, by 
(1 + 0(A/(T -1 )) ^2i(d l r +1 fi) 2 . In contrast with that, the surface integrals of the remaining components of the energy 

momentum tensor arc dominated by J2i ^(<9r<M 2 — (dq^i) 2 ^ ■ This expression can be approximated, in the same limit, 

byO(A/a- 1 )E i (^ +1 /0 2 - 

Thus under the above conditions the electromagnetic radiation undergoes the process of partial isotropisation and 
behaves like a gas with nonisotropic pressure (the tangent pressures Tqq and are negligible) . (That is an interesting 
fact in itself, which in principle opens a way for the explanation of the long standing riddle concerning the apparent 
mass over-abundance of spiral galaxies [10].) 

It was claimed in [1] that all components of the energy-momentum tensor different from T®,T® and T£ must 
vanish due to the isotropy of the discussed cosmological models. Robertson's statement is obviously wrong, but the 
conclusion that T°, Tq and T£ are the only components that matter (in the limit of geometric optics, defined earlier), 
is correct. Therefore the remaining part of the reasoning of [1] can be left intact, modulo another mistake mentioned 
in the introduction (see also below). Let us add here, for the convenience of an inspective reader, that the quantity 
S^ 1 (9' +1 /;) 2 is equal (modulo a constant factor) to the quantity V in the paper of Robertson (see formula (8) in 
[I])- 



A. The proper energy 

The comoving observer will detect the local electromagnetic energy density —T ilv n! 1 n v = —T§. The total energy of 
a wave pulse, 

E(r,t) = - [ dVTl (31) 

JV(r) 

is related to the quantity m by 

The volume V(r) extends over an annulus (ve + rj — A,te + t]), with the areal radial extension a A being of the order 
of the longest wavelength present in the pulse. The energy flux through the comoving coordinate sphere S(r) reads 

4-7T 

8 t E(r,t) = _(^ +1 (.A)) 2 - HE{r,t). (33) 
er 

One finds, assuming (as in the preceding section) that A << a(rE) 1 that the energy term E(r,t) is approximated 
by 4tt f r ^_ A dr^ l {d l + 1 (fi)) 2 /a » (Aa)Zfg 1 (d l + 1 (f t )) 2 / a 2 . It appears that H w l(r 15 /fcm (remember that 
c = G = 1 ) - that value is prohibitively small, so that for any realistic extension A of the radiation pulse the quantity 
AaH « 1. Therefore the second term in (33) is always negligible and the energy flux is well approximated by 

d t E(r,t) = ^(d l r +1 f?J . (34) 

Notice that the factor (<9 r //o) 2 is constant along the trajectory of the outgoing wave pulse. In order to compare this 
with Robertson, let us observe that d t E(r,t) is equal (neglecting the -ff-related term) to 4:nr 2 a 2 U where U is given 
by formula (2.10) in [1]. 
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B. The luminosity distance 



We define (r^,^) and (ro,to) as coordinates of the emitter and the observer, respectively One easily finds, using 
(32) and the constancy of (<9 r /; ) 2 , that 

dtE(r,t)\ to =p^d t E(r,t) \ tE = 

d t E(r,t)\ tE . (35) 



(1 + z) 2 

Notice that dtE(r,t)\t E is the luminosity of the emitted radiation. The apparent luminosity is given by the ratio 
d t E(r,t)\ to /(AttRq). From the definition of the luminosity distance follows now immediately 

D={\ + z)Ro, (36) 

in accordance with the formula (1). 



VI. APPENDIX 

Let / G C°° have a compact support (a, b) such that A = b — a « a. It is easy to see that \\d l f /a\\L 2 ( a r ) << 
\\d l+1 f\\ L 2 (a,r); f° r an y interval (a,r) G (a, b). Therefore the amount m^ a r ^ of radiation within the annulus (a,r) is 
well approximated by |<9 Z+1 /| \ 2 L , , and the pointwise equality of Eq. (30) holds true almost everywhere. 
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